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Abstract 

The exact solution for a system with two-particle annihilation and decoagulation has been studied. 
The spectrum of the Hamiltonian of the system is found. It is shown that the steady state is two-fold 
degenerate. The average number density in each cite (rii(t)) and the equal time two-point functions 
(ni(t)nj(t)) are calculated. Any equal time correlation functions at large times, (nj(oo)nj(oo) ■ ■ ■), 
is also calculated. The relaxation behaviour of the system toward its final state is investigated and 
it is shown that generally it is exponential, as it is expected. For the special symmetric case, the 
relaxation behaviour of the system is a power law. For the asymmetric case, it is shown that the 
profile of deviation from the final values is an exponential function of the position. 



1 Introduction 

In recent years, reaction-diffusion systems have been studied by many people, using different meth- 
ods. Among them are the field theoretic methods, which allow for perturbative approaches to build 
up correlations in low dimensions (lj §]. As mean field techniques can not be used for low dimensional 
systems, people are motivated to study stochastic models in low dimensions, which can be solved exactly. 
Moreover, solving one dimensional systems should in principle be easier. Applying a similarity trans- 
formation on an integrable model, one may construct stochastic models, their integrability may be not 
obvious. Recently, Some people have studied such transformations [3-8]. 

Exact results for some models in a one-dimensional lattice have been obtained, for example in H |9|. 
In These cases, the time evolution of the system is determined by a master equation ]10| . Models with no 
diffusion received less attention in the literature [11-14]: It is said that unless the system has long-range 
reactions ]l2| , |l3| , the time dependence involves exponential relaxation rather than power law behaviour 
typical of the fast diffusion reactions. 

In jlB), a 10-parameter family of stochastic models has been studied. In these models, the fc-point 
equal time correlation functions (riiiij ■ ■ ■ rik) satisfy linear differential equations involving no higher- 
order correlations. These linear equations for the average density (rii) has been solved. But, these set of 
equations may not be solved easily for higher order correlation functions. The spectrum is also partially 
obtained. The model which we address in this article is a special case of that 10-parameter stochastic 
model. 

In this work, we report the exact solution for a system with two-particle annihilation and decoagu- 
lation. This model may be considered as a biased voting model, in the sense that there are two different 
opinions. If the two persons on two adjacent sites have different opinions, they may interact so that their 
opinions become the same. The bias parameter corresponds to the dominance of the left (or right) sight. 



In the absence of bias, this system is equivalent to the zero-temperature Glauber model 16, p.7| . This 
system is related to free fermion system, through a similarity transformation, and hence is solvable. Note 
that the system itself is not a free fermion system and can not be solved by applying only Jordan- Wigner 



transformation. 

When there is right-left symmetry, the average density decays to its final value in the form of power 
law (t~^). But in the general case (biased model) it decays in the form of an exponential. Moreover, 
the profile of the deviation of the average density from its final value is not uniform but exponential in 
terms of the site number. In fact, the parameter representing the right-left asymmetry, in some sense, 
determines the dominance of the right sites over the left sites, or vice versa. 

The spectrum of the Hamiltonian of the system is found. It is shown that the steady state is two-fold 
degenerate. The probability of finding the system in each of these two states is determined by the initial 
average density, and is time-independent. It is shown that at large times, any n-point function is equal 
to the 1-point function, which is position-independent. 

(m(oo)nj(oo) ■ ■ ■ n fe (oo)) = (n^oo)) = -j- ^(n m (0)) (1) 

m 

This is due to the fact that the system has two steady states; either completely full, or completely empty, 
as it will be shown. This means that the mean-field approach does not work and this system is highly 
correlated. 

The scheme of the paper is as follows. In section 2, similarity transformations relating stochastic 
systems to other (stochastic or non-stochastic) systems are investigated. In section 3, a solvable model 
is obtained through a similarity transformation on a free-fermion system. The spectrum of the system is 
also obtained in this section. In section 4, the 1-point function is calculated and its large-time behavior 
is investigated. In section 5, the two-point function and its limiting behavior is obtained. In section 6, 
The null vectors of the Hamiltonian are obtained and from that the steady state of the system is obtained 
in terms of its one-point function at t = 0. Finally, in section 7 we consider the next-to-leading term of 
the one-point function at large times, and from this obtain the way the system relaxes to its final state. 
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2 Similarity transformations as a method for obtaining solvable 
stochastic models 

Here some standard material [2,3,5] is introduced, just to fix notation. The master equation for P(a, t) is 

^P(a, t) = £ [w(t - a)P(r, t) u (a - r)P(a, t)} , (2) 
where lu(t — > tr) is the transition rate from the configuration r to a. Introducing the state vector 

\P(t))=J2P(a,t)\a), (3) 

tr 

where the summation runs over all possible states of the system, one can write the above equation in the 
form 

^\P)=H\P), (4) 

where the matrix elements of 7i are 

(a|W|r) = uj(t — > a), a, 

(a\n\a)=-J2^^r). (5) 

The basis {{cr\} is dual to {|<r)}, that is 

(a\r) = 5„, T . (6) 

The operator is Ti is called a Hamiltonian, and it is not necessarily hermitian. But, it has some properties. 
Conservation of probability, 

£P(M) = 1, (7) 



shows that 



where 



(S\H = 0, (8) 



<s| = £</?!■ (9) 

/3 



So, the sum of each column of H, as a matrix, should be zero. As (5*1 is a left eigenvector of Tt with 
zero eigenvalue, TL has at least one right eigenvector with zero eigenvalue. This state corresponds to the 
steady state distribution of the system and it does not evolve in time. If the zero eigenvalue is degenerate, 
the steady state is not unique. The transition rates are non-negative, so the off-diagonal elements of the 
matrix TL are non-negative. Therefore, if a matrix TL has the following properties, 

(S\H = 0, 

(10) 

(<t\H\t) > 0, 

then it can be considered as the generator of a stochastic process. The real part of the eigenvalues of any 
matrix with the above conditions should be less than or equal zero. 
The dynamics of the state vectors (Q) is given by 

|P(f)> = exp{tH)\P(Q)), (11) 

and the expectation value of an observable O is 

(0)(t) =Y J 0{a)P{a,t) = (S|Ocxp(m)|P(0)). (12) 
tj 

If TL is integrable, one can solve the problem, that is, one can calculate the expectation values. Suppose 
now, that a Hamiltonian is integrable but is not stochastic. There arises a question, whether or not 
there exist exist a similarity transformation which transforms it to a stochastic integrable Hamiltonian. 
Consider an integrable Hamiltonian TL. The similarity transformation 

TL := BUB- 1 (13) 

leaves its eigenvalues invariant. Consider a special case: The system consists of a one dimensional lattice, 
with nearest-neighbor interaction, 

L 

H = ^W M+ i. (14) 

i=l 

Suppose, also, that the system is translation-invariant: 

(15) 



and we are using periodic boundary conditions. A simple class of similarity transformations is then 

B = r 1( g)r 2 (g)---®r L . (16) 

The simplest case is when all IYs are the same. In this case, if one can find T such that 

H = r®THT- 1 & (17) 

is stochastic, then TL defined through ( |l3| ) would be stochastic. A more general class of similarity trans- 
formations is obtained through 

T t :=T{g)\ (18) 

where g should have the property 

[g®g,H]=Q. (19) 

In this case, one obtains 

H = {T(g>T g) H(T(g>T g)- 1 . (20) 
Define (s\ to be the sum of all bra-states corresponding to a single site. We then have 

(S\ = (s\ <g> • • • ® (s\ . (21) 

L 

For H to be stochastic, its off-diagonal elements should be non-negative, and we must have 

(s\ <g> (s\H = 0. (22) 

This shows that 

H®(£| :=(s|r® (s\Tg, (23) 

should be an eigenvector with zero eigenvalue of H, that is, H should have a decomposable left eigenvector. 
So, in order that this prescription of constructing integrable stochastic model works, one must begin with 
a Hamiltonian H, the left eigenvector with zero eigenvalue of which is decomposable. The real part of 
all other eigenvalues of H should, of course, be non-positive. 
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3 A one— parameter solvable system on the basis of a free— fermion 
system 



Consider the Hamiltonian 



sj +1 s t -ni{l-n i+ x)\ 



i=l 

1 — r) . 



U+i s t -n i+ i(l-m)] 



+H s i+i s i -n t n i+1 }}, 



(24) 



where s + , s , and n are 



1 















: 


')■ 










0/ 






0/ 



(25) 



and the subscript i represents the site, on which the operator acts. This Hamiltonian describes the 
following processes 

1 + 7? 



A0 -> Via 
A® 

AA —> 00 



with the rate- 



with the rate 



2 

1-T) 



with the rateA. 



(26) 



This model has been recently studied. In the case A = 0, the above model describes an asymmetric 

exclusion process. For A = 1, the Hamiltonian is bilinear in terms of creation s + and annihilation s~ 

operators. This problem has been solved via a Jordan- Wigner transformation p8[ |l9| |. In the notation 

of the previous section the matrix form of H is 

/—A 0\ 

_1±2 i-Ii 



H := 



1+3 

2 



l-»7 







(27) 





V A oJ 

This matrix has two eigenvalues, and -1, both of them are two-folded degenerate. One of the zero 
left eigenvectors can be decomposed into a tensor product. Doing the above mentioned procedure, this 
Hamiltonian can be transformed to another stochastic one. For this case, One can show that the matrix 



g is the identity matrix, and the similarity transformation for all sites become the same. This has been 
done in ||. 

One of the left eigenvectors corresponding to the eigenvalue -1 has also the desired property. To use 
the prescription described in the previous section to construct a stochastic Hamiltonian, we define a new 
Hamiltonian, 

H':=-H-l, (28) 

and apply the similarity transformation on this new Hamiltonian. One of the zero left eigenvectors of H' 
is 



(1 0) = (1 0)®(1 0). 



The similarity transformation should map (s\ <X> (s| to this vector: 



(29) 



(1 l)r®(l l)Tg = a(l 0)0(1 0) 



(30) 



(31) 



So, 

(i i)r = ai/(i o) 
(i i)r fl = 2(i o). 

Scaling T and g does not alter the Hamiltonian H . So we can remove a and v by scaling the matrices 
r and g. Then the above relation gives some constraints on the elements of T and g. the condition of 
positivity of rates, fixes g and T: 
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1 -1 

1 1 
1 



(32) 
(33) 



-1 

The two site Hamiltonian, then, takes the following form 

0-100 
0-10 



VO i±2 i=2 0/ 



H = 



(34) 
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and TL is 



i=l 

1 + V r + - t-, \n 

+ — — [sjn t+1 + s i (1 - n i+ i)J 

- [nj(l - n i+ i) - (1 - Uj) rii+i]}. 

This Hamiltonian describes the following processes 

A®-> AA 
A0^00 
®A^ AA 
0A-> 00 



1 




'/ 




2 




1 


+ 


'7 




2 




1 


+ 


'7 




2 




1 




V 



(35) 



(36) 

The Hamiltonian ( |35| ) is not quadratic in s + and s — . So, one can not map this Hamiltonian to a free 
fermion system, using a Jordan-Wigner transformation. But the Hamiltonaian TL is integrable and can 
be mapped to a free fermion system by a Jordan- Wigner transformation. Consider the following Jordan- 
Wigner transformation (fijl, [l9| 



,t .- 



Qj-iSj 

, ■ Qi-i't 
Qj ■ ri(-4)- 

It can be easily shown that the number operator at each site is, in terms of new generators, 

1 + s? 



(37) 



ni := 



(38) 



Using this transformation, The Hamiltonian TL takes the following form 

nj sr l" 1 ~^ t , 1 + 7 ? t t 



(39) 



and gJ fulfill the fcrmionic anti-commutation relations 



{fli,aj} = {a], a]} = 



{ai,a|} = Si- 



(40) 



Note that it is in the limit L — *■ oo that the Jordan- Wigner transformation we are using, works. Otherwise, 
there are some boundary terms in (|3^) as well. So, all the results we obtain hereafter, are valid only in 
this limit. Now, introducing the Fourier transformation 



1 v-^ , t , — 2ni j k. 



and substituting it in (MOT), it is seen that 



(41) 



{b k ,b l } = {blb]} = 



{b k M} = S. 



(42) 



As a result, the Hamiltonian H takes the form 



n 



E 
E 



1 — 7] lixi k 1 + rj — 2ni k 
— exp( — ) + — exp( — — ) 



£kb\b k - ism(^^)b^ k b k 



blb k + 6- fe 6 fc exp(- 



blb k 



(43) 



where 



27rfc 27rfc 
£fe := -1 + cos(-^— ) - z 77sm(-j— ) 



AO 



One can now, easily obtain the time dependence of b k and bj., using ([42|) and — ^- = [O, H] 



(44) 



hfc(t) = ^>fc(0)e £ 



ftt ( t) = e-«»*{&t (0 ) - JC ot(— )[ e («*+ e -*)* - l]&_ fc (0)} 



(45) 



Now we return to our problem: determining the expectation values of a system evolving with the Hamil- 
tonian H. The expectation value of a quantity O is 



(0)(t) ={S\Oexp(tH)\P(0)) 

= (S\exp(-tH)Oexp(Ht)\P(0)}. 

Substituting H — —BHB^ 1 — LI, where 1 stands for the identity matrix, yields 



(46) 



(0){t) = (tt\6(-i)B^\P{0)) 



(47) 



where 



6 —B^OB, 



tH/r\„—VH 



0(-t) :=e tn Oe 



and 



(fi|:=(l 0)0(1 0)O---(l 0). 



The main expectation values of interest are the correlation functions of n/s. To determine these, 

r- 1 nr = i(i-s+- s -) 
(r 9 )- 1 nr g = i(i + s + + s -) 



So 



B- 1 n t B=-[l-(-iy( S + + S -)]. 



Now, we want to calculate the expectation value of O where 



We have 



O ■■= n im ■ ■ -n^rii^ h{ii{- ■ ■ (i r , 



° = ^ [i - (*l + s- m )]■■■[!- (-ir (4 + ^ )] • 



Using the Jordan- Wigner transformation, one arrives at 



^ = ^[ 1 -(- 1 ) im ft m -i(«L+«i 



It is easy to check that (fl\Qi = (— l) l (^|- So in calculating (£>), one can use O' instead of O: 



Instead of C, It is enough to set O" in the expectation value of O, where 



O" :=^(1 + «L) - ■(! + <)■ 



To prove this, one should use (fi|7i = — L(fi| and (f2|aj(0) = 0. 
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4 The one— point function 

As the first example, consider the one-point function (n m (t)}: 

(n m (t)) = l<n|[i + a t l (-*)]B- 1 |P(Q))- (58) 

Using the Fourier transformation (pi]), the time dependence of b' k (p5[), and remembering (f2|6fc(0) = 0, 
we obtain 

(n m (t)) = ^ + ^=^e"^(fi|^(0)S- 1 |P(0)) e ^ t (59) 



1 1 

T ' 

k 

Now, we use the inverse Fourier-, and Jordan- Wigner-transformations, and arrive at 



(%(*)) = \ + ^ E e? ^^ (S\Bsp- l \P{Q))e*»\-iy-\ (60) 
This can be written in a simpler form, using 

. . 2n, — 1 + s7 — sj" 

Ba+B- 1 = {-I) 1 - 1 — ^ 3 - (61) 

and 

(s|(2n- 1) = (s\(s~~ - s+). (62) 

One then arrives at 

W*)) = 5 + ^E f ^^<S|(2n,(0) - l)l^(0)) e £fct . (63) 
Using (S|P(0)) = J2a P ( CT ' °) = !> one arrives at 

(%(i)) = E A ™i( t )("i(°)>. ( 64 ) 

where 

A mj (i):= z E e ^^ e£fc ■ ( 65 ) 
Now, consider the limit t — > oo. In this limit, the only contribution in the above summation comes from 
the term k = 0. So, 

lim<n ro (*)> = yEv^°)>' (66) 

3 
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which shows that in the limit t -> oo, the expectation value of the number of particles in any site tends 
to the average of the initial value of this quantity. In the last section, we will find the next leading term 
of (rii(t)), for large times. 

Now, we want to calculate the expectation value of the number of particles in the site j in the limit 
L — > oo . First, we calculate A m j in this limit. To do so, we define z :— exp(i^J^). We then (in this limit) 
arrive at 

A^W^e-'^^z^exp^i^z + ii^z- 1 )]. (67) 
Changing the variable z to w := zwjij, the above integral takes the form 



or, using the change of variable w :— e l8 , 



k mj {t) = e-t 2 ^ e *W-m)o+*Vi=^co.^ (69) 

The above integral is an integral representation of the modified Bessel function: 

W*)=(l^) 2 lrn- 1 (tVT^)e- t . (70) 



(n m (t)}, m the limit L — > oo, is then 



Mt))^j2{hr) ^-^v/wk'mo)). (7i) 

5 The two— point function 

The other quantity which we want to calculate is (n m (t)ni(t)) . Without loss of generality, one may 
assume (m > I). To calculate this, we use (|^), which gives 

(n m (t)m(t)) = \(n\[l + al(-ty[l + a}{-t)]B- 1 \P(0)) 

= -| + l[{n m (t)) + (n,(*))] + l{n\al(-t)al(-t)B-i\P(0)). (72) 
The main thing is to calculate the last term. To do this, we first use the Fourier transformation of aj's, 
(QlaK-ilaU-fyB-'im) = tE 6 "' (^(-^H")^ 1 !^)), (73) 



L 

k,p 
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and then substitute the time dependence of bts. 

(fi|ot l (-t) a t(-t))B- 1 |P(0)) = ±X>^^^+C e *+^'<fi|&;(0) 

k,p 

[ftt(O) + icot(^) (l - e-^-^) 6_ p (0)] B^IPCO)). (74) 
Now we use inverse Fourier transformation for the term. The other term is easily summed. We 



arrive at, 



r,s 

fc 

or, 

(n\al(-t)al(-t)B^\P(0)) = A mr (t)A Ja (t)((2n r -l)(2n.-l)) sgn(r-a) 

+ ± £ e i cot( 7rfc } ^ _ e(£fe+£ _ fe )^ 5 (?6) 
k 

where we have used the definition of A^, and (• ■ -)o means the expectation value at the initial time. 
Adding all terms in (jT^) together, one arrives at 

(n m (t)m{t)) = -j((n m (i)) + (m(t))) +^ A mr (t)A ls (t)((n r - i)(n a - ^)) sgn(r - s) 

+ ^ E el (l " e^ + ^) . (77) 

fe 

The last term is independent of initial conditions, So we can calculate it for a special case, e.g. \P(0)} = 
|0). Then, the final result is 

(n m (t)m{t)) = -((n m (t)) + (nj(t)» +E A ror (t)A is (t)sgn(r - s)(n r n s - ^y^)o (78) 

For large times, it is seen that 

lim<nm(t)ni(t)) = <n(oo)) (79) 

f — »oo 

6 Null vectors of the Hamiltonian, the steady state of the sys- 
tem, and the n— point function 

Now we want to study the null eigenvectors of the Hamiltonian H. It is easy to see that the Hamiltonian 

13 



( pig ) has at least two null eigenvectors, which means that the steady state is not unique. These states are 
one in which all sites are occupied, and one in which no site is occupied. One can check this easily by 
acting the Hamiltonian (p5|) on these states. It was shown that the Hamiltonian H may be written as 



n 



E 

k 

E 

fe>o 



£kb\b k - isin(— — )b^ k b k 



e k b\b k + e_ fe 6L fc 6-fe - 2 i sm(^^)b- k b k 



(ob\b Q . 



(80) 



This Hamiltonian is obviously block diagonal. In each four dimensional block, one can choose a basis 
{|0), &|.|0), 6_ fc |0), b k b^_ k \G)}. The eigenvalues of this four dimensional block are 0, e k , e_ fc , e k + 6- k , or 
A/fcCfc + Af-kC-kt where TV's are zero or one. The eigenvalues of the Hamiltonian are, therefore, 



E{N} = Y J M k e k . 



(81) 



From this, one can obtain the eigenvalues of Ti as 



Here we have used 



*W = -S(JV* e fc + 1 ) 

k 

= ]T(-A/" fe + l)e fc . 

k 

$> + e*) = 0. 



(82) 



(83) 



Now, it is easy to sec that E is zero iff 1 — J\f k = 0, Vfc ^ 0. This shows that the null eigenvector is 

two-fold degenerate. As the final state is two-fold degenerate, and it is known that the totally full state 
/1\ /1\ fl\ f0\ f0\ /(T 

:= I ® (g) • • • ) and totally empty state (|0) := I <g> <g> • • • ) ) are null 

W W W V 1 / W V 

eigenvectors of the system, we have 



where a + (3 = 1. Using (|66|), and 



it is seen that 



|P(oo))=a|0)+/?|fi>, 

(S\n t \0) = 0, 
(S\rn\Sl) = 1, 

^(nH)4E(%(0)) 



(84) 



Po 



(85) 



(86) 
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and 

a = 1 - po. (87) 

From this, one obtains 

|P(oo)) = [l-p ]|0) + p |O). (88) 
Using this, it is easy to find all m-point functions in the limit t — > oo. We have 

(rii m •••n il )(oo) = (5|n im • • •n il |P(cx))) 
= Po 

= ("M) 

(n(oc)) m (89) 
This clearly shows that the mean-field approximation does not work here. 

7 Relaxation of the system toward its steady state 

It was shown that in the limit t — > oo, the expectation of the number of particles in any site tends to the 
average of the initial value of this quantity. Now, we want to study the behaviour of the system at large 
times. Starting from (fFl|), and representing (rij(0)) by its Fourier transform, we have 

(n ro (t)) = e- t jf W ^^^^y m _ j (t^/T^)e- i ^[f(u) + 27rn5(u)}, (90) 

where [f(u) + 2irn6(u)] is the Fourier transform of (n^O)), and n is the average density. We have extracted 
this part of the Fourier transform, so that the remaining is a smooth function of u. Then, f(u) denotes 
the Fourier transform of the deviation (rij(0)) — n. The summation on j is easily done, using 

^2x n l n (y) = e {y/2)(x+1/x) . (91) 



So. one arrives at 



(n m (t)) = fi + pe^fi^e^, (92) 
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where e(u) is the same as eu with k = Lu/(2n). The above integral is simplified for large times, using 
the steepest descent method. Using the change of variable z :~ e lu , the integral becomes 

(n m (t)) =n+l ^- e ^- 2+ ^ z+ ^ z ~^ 2 f(z)z- m , (93) 
J ziriz 

where the integration contour is the unit circle. The multiplier of t in the exponent is stationary at 

* - M < m > 

and 

As the real part of this multiplier is larger at z = z\ , the integral gets its main contribution from this 
point. (This point is not on the integration contour. But, assuming / to be analytic, one deforms the 
integration contour so that it passes from zi, and then uses the steepest descent method.) We arrive at 

<n m (i))-n~i=(i^jj e «(V^-i). (96) 

The effect of the Fourier transform /, and the second derivative of the multiplier of t in the exponent is 
a multiplier independent of m and t. Two general features, independent of the initial condition, are seen 
from the above relation. First, the decay to the final state is not in the form of a power law, but in the 
form of an exponential. It becomes a power law only in the symmetric case r\ = 0. Second, if r\ > 0, the 
expectation at the rightmost sites tends rapidly to its final value. That is, the profile of the deviation 
from the final value is decreasing with respect to m. This is so, since in this case the two-site reaction 
is favorable to the state where the left site changes so that it becomes identical to the right site. This 
means that cases where the right-site changes is less probable than cases where the left site changes. So, 
the right site arrives earlier to its final state. This expression seems to be unbounded for either m —* oo 
orm^ — oo. For any fixed t, this is true. But it simply means that in order that this term represents 
the leading term for some m, t must be greater than some T, which does depend on m. 
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